Strain-controlled correlation effects in self-assembled quantum dot stacks 
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We show that elastic interactions of an array of self- assembled quantum dots in a parent material 
matrix are markedly distinct from the elastic field created by a single point defect, and can explain 
the observed abrupt correlation-anticorrelation transition in semiconductor quantum dot stacks. 
Finite volume effects of the quantum dots are shown to lead to sharper transitions. Our analysis 
also predicts the inclination angle under which the alignment in successive quantum dot layers occurs 
in dependence on the material anisotropy. 

PACS numbers: 68.65.Hb,81.16.Dn,81.16.Rf 
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The self- assembled growth of semiconductor quantum 
dots (QDs) is widely used as an efficient tool for fabri- 
cating nanostructures with important potential applica- 
tions in opto- and nanoelectronics Q, 0. The process 
of self-assembly is based on the formation of coherent 
localized structures confined in three spatial directions, 
i.e. QDs, on a flat, uniform two-dimensional (2D) wet- 
ting layer in highly strained materials, using the Stranski- 
Krastanov growth mode. Multilayer systems in the form 
of vertical stacks are used to improve the lateral ordering 
of QD arrays, and minimize the size and shape fluctu- 
ations of the QDs. By tuning the spacer thickness a 
certain desired lateral spacing of QDs can be obtained, 
and one can switch between several types of stacking 
modes Q Q Theoretically, several approaches have 
been made to understand the new features of QD stacks. 
Elasticity theory has proven to be very useful, not only 
for understanding the growth of QD systems consisting 
of one layer 7], in particular when combined with ki- 
netic Monte Carlo simulations 0, 0, [HI, but also for 
givin g in sight into the growth of stacked QD systems 
E 13 Hi El H 113 0ur understanding of the verti- 
cal and lateral correlation and anticorrelation properties 
of such stacked arrays of QDs is, however, still incom- 
plete. Recently, clear experimental evidence of an abrupt 
transition between vertically aligned (correlated) and an- 
tialigned (anticorrelated) QDs in the InGaAs/GaAs ma- 
terial system in dependence upon the spacer layer thick- 
ness has been presented Q , but this transition could not 
be explained by the anisotropic strain field of a single QD 
0, E 5 ari d the angle against the vertical under which the 
alignment in successive QD layers is observed seriously 
disagrees with the theoretically predicted tilt angle for 
the minima of the strain |5|. 

In this report we present a three-dimensional model 
of the strain field of a stacked array of quantum dots, 
taking fully into account the anisotropy of the material 
system and the 3D shape of the QDs. This is an ex- 
tension of the model proposed by Shchukin et al. [T^ . 
While other theoretical descriptions 0,0 using a similar 



Green's tensor approach [T]| have employed strictly pe- 
riodic arrays and point-like or stripe-like approximations 
of the QDs, we have implemented arbitrary QD shapes 
and configurations. Within our approach we are able to 
explain both the experimentally investigated tilt angles 
and positions in an anticorrelated stacked array of QDs 
in the InGaAs/GaAs material system |fj and describe 
the transition from correlated to anticorrelated growth. 

We model the multi-sheet array of QDs by inclu- 
sions of material A in a matrix of material B. Starting 
from elasticity theory for anisotropic materials 0, ^| , 
we solve the equilibrium equation Vj (A^Vfci^(r)) = 

Vj(cr^i?(r)) with the elastic moduli Xijkh the elastic dis- 
placement field ui, the characteristic shape function i?(r) 
of the QDs ($(r) = 1 if r is inside an inclusion and zero 
elsewhere), and the stress tensor aff which is connected 

,(0) Kt: „(o) 



to the stress-free strain tensor e\y by a\j J = \jkl^ki 
Here Einstein's summation convention is used. Since the 
volume of the inclusions is assumed to be small compared 
to the parent material, the homogeneous moduli approx- 
imation is used. The problem is solved under stress-free 
boundary conditions on the surface. 

In general, the elastic energy of interacting inclusions 
is given by 
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Introducing the in-plane coordinates rn and the vertical 
coordinate z, one can express the elastic interaction en- 
ergy per unit surface area in the growth plane at z = zo 
in terms of the static Green's tensor as [l2j 

E(v f z = z ) = h^ f dVy f ° dz'#{r' r ,z , )W{r n -r\;zo-z') 

J J — oo 

(2) 

where the integral is over the buried QDs with z r < zo, 
is the height of one single inclusion on the surface, 
and the interaction term 
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is related to the Fourier transform of the static Green's 
tensor On, which can be computed numerically for elas- 
tically anisotropic cubic crystals, for details see [13, 



tically 



Z (c) 


7(b) 


Z ,a, 

z' 



(a) 



(b) 



(c) 



IH1 











"D 




— 


o- 





FIG. 1: Schematic representation of the dependence of the 
surface elastic energy upon the spacer thickness. The cross- 
sections of samples with different spacer thickness d (a)-(c) 
are combined into a single cross-section (d). 
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FIG. 2: Elastic strain field of one single point-like inclusion 
(a) and a periodic array of point defects (b)-(d) in a GaAs 
matrix, (a) and (c) show the minima (black) and maxima 
(grey) of the elastic energy E per unit area in the [110] di- 
rection of the relaxed surface vs. spacer thickness d. The 
density plot (b) depicts E in the [001]- [010] surface plane at 
d = 25. (d) E for d = 10 and d = 25. The positions of the 
QDs are indicated by triangles (a),(c),(d) and by boxes (b), 
respectively. All lengths are in units of the lattice constant of 
the parent material aGa As = 0.565 nm (monolayer width). 

We have computed the elastic energy in the relaxed 
surface at z = zo for a given array of buried QDs at 
z = z' in dependence on the spacer thickness d — zq — z' . 
The results are combined into a single cross-sectional plot 
as schematically illustrated in Fig. ^ Periodic boundary 
conditions in the lateral direction are used. Fig. Et a ) 
shows the minima (black) and the maxima (grey) of the 
elastic interaction energy E in the relaxed surface for a 
single buried point-like QD. The minima and maxima oc- 
cur in the [110] direction at angle a = 25° and a = 60° 
against the vertical, and correspond to maximum tensile 
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FIG. 3: Elastic strain field of a periodic array of pyramids 
in a GaAs matrix, (a) density plot of the elastic energy per 
unit area E in the [001]- [010] surface plane at d = 35 (inset: 
3D view) (b) minima (black) and maxima (grey) of E in the 
[110] direction of the relaxed surface vs. spacer thickness d. 
(c) E for d = 15 and d = 35. The positions and sizes of the 
QDs are indicated by boxes (a) and dotted triangles (b),(c). 



strain (E < 0, black) and maximum compressive strain 
(E > 0, grey), respectively. Regions of tensile strain give 
rise to an attractive potential and act as preferred nu- 
cleation centers for QDs, whereas regions of compressive 
strain form repulsive potentials pj| . 

As can be seen, our model is in good agreement with 
previous theoretically predicted values of the angle a in 
which the minima of the elastic energy in the [110]- [001]- 
plane occur |4j. Computing the inclination angle a for 
other material systems, we find that a increases with 
increasing anisotropy ratio A = 2c^/(cn — C12) of the 
parent material, where en, C12, C44 are the elastic moduli 
in Voigt notation. We obtain a = 19° for Si (A = 1.56), 
a = 25° for GaAs (A = 1.83), and a = 33° for ZnSe 
(A = 2.04), in agreement with However, these results 
contradict the experimental observation of much larger 
angles p. 

The situation changes drastically if we consider an ar- 
ray of QDs (Fig. EJb)-(d)). We use a quadratic array of 
4x4 QDs with a lateral spacing of I = 25 lattice sites 
(l.s.). In a neighborhood of approximately 1/2 around 
each QD the elastic properties are governed only by this 
single point defect. Outside this region the strain fields 
are overlapping, which leads to new effects (Fig. Etc)). 
Starting with a spacer thickness of d ~ 14 l.s., the inclina- 
tion angle a increases more and more until at d = 17 l.s. 
the two minima induced by two neighboring QDs meet 
and form a flat double minimum exactly in the middle 
between these QDs. With further increasing spacer thick- 
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ness the position of the double minimum remains stable 
and unchanged for a considerable range of d-values. Fi- 
nally, at d = 31 l.s the minimum starts moving back to a 
position vertically above the QDs. 

The transition of the minima from a position close to 
vertically above the QDs to in between the QDs is also 
visible in the elastic energy profiles at d = 10 l.s. and 
d = 25 l.s. in Fig. Efd). This indicates a transition 
from correlated to anticorrelated growth. From Fig. [2(c) 
we infer an inclination angle for alignment of QDs of 
about a = 45° for GaAs which is in reasonable agreement 
with the experimentally observed value of 50° 6]. The 
small remaining discrepancy might result from several 
subtle effects. First, in our model the material of the 
QDs is completely neglected. Second, stacking of the 
QDs might also be influenced by modulations in chemical 
composition or be due to morphological changes in the 
spacer layer. 

For different materials, again the inclination angle a 
increases with increasing anisotropy, as in the case of 
one single QD: a = 43° for Si, a = 45° for GaAs, and 
a = 47° for ZnSe. 

We have also studied the influence of the shape and 
volume of the QDs. In Fig. 03 we consider a 2 x 2 ar- 
ray of pyramids, each with a baselength of 20 x 20 l.s. 
and a height of 5 l.s positioned with a lateral distance of 
I = 40 l.s. Compared with an array of defects, there are 
significant differences. The inclination angle a remains 
unchanged. But the finite volume of the QDs results in a 
large range of spacer thickness d where the minimum en- 
ergy is directly located vertically above the buried struc- 
tures (Fig. El (b)). At a spacer thickness of about four 
times the height of the pyramids the energy minimum 
splits and moves towards the anticorrelated positions be- 
tween the QDs. Due to this finite size effect, the tran- 
sition from the correlated to the anticorrelated regime 
happens much more abruptly than in an array of point 



defects (Fig. Etc)). Such an abrupt transition is indeed 
observed experimentally [(j. Also the double minimum 
of the energy profile is replaced by a single minimum 
(Fig. 0(c)) • Note that the transition occurs at larger 
spacer thickness d « 32 for the array of QDs with a 
larger lateral distance of I = 40 l.s., and this ratio should 
scale for arrays with even larger lateral distance at fixed 
angle a. This is consistent with the experiment p where 
the transition was observed at a spacer thickness of 150 
monolayers for an in-plane nearest-neighbor distance of 
80-100 nm. 

In conclusion, we have shown that elastic interactions 
of an array of QDs in a parent material matrix explain the 
observed correlation-anticorrelation transition, in con- 
trast to calculations which take into account only a single 
point defect. The observed material trend of increasing 
inclination angles of the vertical stacking mode for larger 
anisotropy is also corroborated. The abrupt transition 
has been attributed to the finite volume of the QDs, 
and the reclining angle under which the anticorrelated 
QDs align is in good quantitative agreement with the ex- 
periment on self-organized InGaAs/GaAs quantum dot 
stacks. 

The algorithm introduced in this paper allows for fast 
computation of the elastic interaction energy of nano- 
structures with arbitrary shape and spatial arrangement 
in three dimensions, and the only parameters are the 
elastic moduli of the considered anisotropic cubic parent 
material. Thus it is appropriate for implementation into 
kinetic Monte Carlo simulations [l(J, [l3[ which might be 
used for a detailed investigation of the growth kinetics of 
self-assembled QD stacks in future research. 
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